Generalized Limits for Single-Parameter Quantum Estimation 
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We develop generalized bounds for quantum single-parameter estimation problems for which 
the coupling to the parameter is described by intrinsic multi-system interactions. For a Hamil- 
tonian with fc-system parameter-sensitive terms, the quantum limit scales as 1/N where N is the 
number of systems. These quantum limits remain valid when the Hamiltonian is augmented by 
any parameter-independent interaction among the systems and when adaptive measurements via 
parameter-independent coupling to ancillas are allowed. 
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Many problems that lie at the interface between 
physics and information science can be addressed using 
techniques from parameter estimation theory. Precision 
metrology, timekeeping, and communication offer promi- 
nent examples; the parameter of interest might be the 
strength of an external field, the evolved phase of a clock, 
or a communication symbol. Fundamentally, single- 
parameter estimation is a quantum-mechanical problem: 
one must infer the value of a coupling constant 7 in the 
Hamiltonian _ff 7 = fvyho of a probe system by observ- 
ing the evolution of the probe due to H 1 0, H S 0, II] • 
We take 7 to have units of frequency, thus making ho a 
dimensionless coupling Hamiltonian. 

Quantum mechanics places limits on the precision with 
which 7 can be determined. It is now well established, 
via the quantum Cramer-Rao bound [l], @, H, 5, that 
the optimal uncertainty in any single-parameter quan- 
tum estimation procedure is ^7 ~ 1/y^t Afep, where v 
is the number of independent probes used, t is the evo- 
lution time of each probe, and Aho is the standard de- 
viation (uncertainty) of ho 0, 13| ■ The 1 / \fv dependence 
is the standard statistical improvement with number of 
trials; generally, for nonGaussian statistics, the sensitiv- 
ity 1/y/vtAho can only be attained asymptotically for a 
large number of trials. Besides increasing the number of 
trials, there are two other obvious ways to improve the 
sensitivity: (i) the probe can be allowed to evolve un- 
der Hry for a longer time t; (ii) the quantum state of the 
probe can be chosen to maximize the deviation Aho- I n 
all practical settings, decoherence or other processes limit 
the useful interaction time; in addition, temporal fluctu- 
ations in 7 often prevent the evolution time from being 
arbitrarily extended. For a given parameter estimation 
problem, ho is fixed, as is its maximum deviation. 

Quantum mechanics does, however, provide another 
opportunity: gathering N probe systems into a single 
probe, which is prepared in an appropriate entangled 
state; if Aho for the entangled state increases faster than 
y/~N, the sensitivity improves, provided there is still a suf- 
ficient number of probes to reach the asymptotic regime 
in number of trials. This Letter focuses on how S'y scales 
with N, the number of systems in a probe. Thus we 



work throughout with bounds on the sensitivity of a sin- 
gle probe, remembering that the bounds can only be 
achieved by averaging over many probes, but preferring 
not to muddy the discussion by carrying along the 1/ 1 yfv 
dependence on the number of probes. 

For the iV systems in a probe, it has been traditional 
to consider Hamiltonians of the form 



H 1 = hjho 



N 



(1) 



where the /ij's are single-system dimensionless coupling 
Hamiltonians, assumed to be identical. Restriction to 
Hamiltonians that are separable and invariant under par- 
ticle exchange, as in Eq. ([T]), is physically motivated: in 
metrology it is generally desirable to make the coupling to 
the parameter homogeneous, and multi-body effects are 
typically undesirable because they are less well character- 
ized. In atomic clocks, for instance, much experimental 
effort is directed toward achieving a Hamiltonian of the 
form ([1]). In many cases, even the measurements per- 
formed on the probe are unable to distinguish between 
individual constituents. 

To determine how the optimal parameter uncertainty 
scales with N, one maximizes the deviation Aho over 
joint states of the probe systems. If entanglement is not 
allowed, the probe systems can be regarded themselves 
as independent probes; in this case, Aho scales as %/iV, 
producing the so-called shot-noise limit found in precision 
magnetometry, gravimetry, and timekeeping [6j. When 
entanglement is allowed, however, one can choose the 
initial probe state to be the "cat state," 



|Amj ■ ■ ■ j Am) + |A m , . . . , A m )j 



(2) 



where for system j, \Xm) (|A m )) is the eigenstate of hj 
with maximum (minimum) eigenvalue Am (A m ). This 
yields a deviation Aho = N(\m — A m )/2 that scales lin- 
early in N [3, |5| , a scaling known as the Heisenberg limit. 
Evolution under H~ t for time t introduces a relative phase 
e w) into the cat state, with 4>{t) = jtN(XM — A m ), and 
leaves Aho unchanged; the Heisenberg limit can be at- 
tained (asymptotically for many trials) by measuring on 
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each probe system an observable two of whose eigenvec- 
tors are |±) = (|A M > ± |A m ))/^. 

The Heisenberg limit is not general since there are 
physical systems of interest for parameter estimation 
where a coupling Hamiltonian of the form ([T]) is overly 
restrictive. In particular, some condensed and even 
quantum-optical systems exhibit nonlinear collective ef- 
fects due to multi-body or tensor-field interactions 0, H, 
9] . In such systems, multi-body terms in the Hamiltonian 
can couple to metrologically relevant parameters. In this 
Letter we generalize single-parameter quantum estima- 
tion to intrinsic many-body interactions, which surpass 
the conventional Heisenberg limit. Our work is lar gely 
inspired by a recent paper of Roy and Braunstein [lOj, 
which claimed an exponential scaling for a collection of 
N qubits with a particular Hamiltonian. We argue below 
that the exponential scaling is unphysical. 

We turn now to showing that Hamiltonians with in- 
trinsic fc-body terms generate a family of parameter es- 
timation problems, characterized by k, where the quan- 
tum limit scales as 1/N k . For this purpose, we consider 
Hamiltonians of the form 



H y (t) = hyho + H(t), h = J2 h hl,j k 



(3) 



where ho is the dimensionless Hamiltonian that describes 
coupling to the parameter. The auxiliary Hamiltonian 
H(t) is discussed below. In ho, k denotes the degree 
of multi-body coupling, with the sum running over all 
subsets of k systems. We could also include couplings 
of different degrees up to a maximum degree, but since 
the maximum degree dominates the sensitivity scaling, 
we stick with a single degree k in the following. We as- 
sume that the fc-body coupling h^ is symmetric under 
exchange of probe systems. Moreover, we assume that k 
and h^ are independent of the number of probe systems. 
We make this latter assumption, that ho is an intensive 
property of the probe, because we want to consider a par- 
ticular kind of coupling to the parameter which remains 
unchanged as N changes. For real physical systems, the 
symmetry and intensive assumptions will hold only ap- 
proximately and only over some range of values of N. 

The auxiliary Hamiltonian H(t) includes all 
parameter-independent contributions to _ff 7 . For 
example, it includes the free Hamiltonians of the probe 
systems and any parameter-independent interactions 
among them. In addition, we can introduce an un- 
determined number of ancillas and let H include the 
couplings of the ancillas to the probe systems and 
any couplings among the ancillas. Measurements on 
the ancillas can be included as part of an overall final 
measurement on the probe-ancilla system; since the 
Cramer-Rao bound that underlies our analysis holds for 
all possible measurements and ways of estimating 7 from 
the measurement results, the bounds we derive hold for 
arbitrary measurements on the ancillas. This conclusion 



applies even to measurements on the ancillas that are 
carried out during the evolution time and whose results 
are used to condition measurements on other ancillas or 
to control the coupling of other ancillas to the probe. 
By the principle of deferred measurement ll|, which 



is illustrated in Fig. [TJ all such measurements can be 
shuffled to the end of the evolution time by making 
appropriate adjustments to H. 

Now let po be the initial state of the probe and any 
ancillas. After a time t, the state evolves to p 7 (i) = 
Uj(t)poU^(t), where the evolution operator is generated 
by the Hamiltonian ([3]): 



in — - ». >./ n 



= H^(t)U 7 (t) 



(4) 



At time t, measurements are made on the probe and 
ancillas, the results of which are used to make an es- 
timate 7 cs t of the parameter. The appropriate statistical 
measure of the estimate's precision is the units-corrected 
mean-square deviation of 7 cs t from 7 0, 0] : 



<5 2 7 
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M(7est)/d7l 



(5) 



Here and below expectation values are evaluated with 
respect to Py(t). 

The quantum Cramer-Rao bound states that [J 0,0,01 
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X 7 (t)=tr(p 7 (t)£ 7 (t)) ={£*(t)) , (6) 



where X 7 (£) is the quantum Fisher information. The Her- 
mitian operator £ 7 (i), called the symmetric logarithmic 
derivative, is defined (implicitly) by 



d p 

(£ 7/ o 7 + p 7 £ 7 ) = = — i\K 1 ,p 1 ] 



Here 



97 

. dU 7 (t) 

d-f 



(7) 



(8) 



is the Hermitian generator of displacements in 7. If there 
is no auxiliary Hamiltonian, Kj(t) = tho. 

For pure states, differentiating p 7 = p 7 shows that 

£ 7 (t) = = -2i[K 7 (t), Pj (t)] . (9) 

Then the Fisher information reduces to a multiple of the 
variance of K^(t): 

J 7 (i) = ±((K*(t)) - (K,(t)f) = 4A 2 K,(t) . (10) 

For mixed states, the variance provides an upper bound 
on the Fisher information, instead of equality 0]. 

We define the operator semi-norm \\H\\ of a Hermitian 
operator H as \\H\\ = Mh — win, where Mh (ttih) is the 
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ancilla 



FIG. 1: Quantum-circuit diagrams illustrating the principle of deferred measurement. In the circuit on the left, a measurement 
M on the lower ancilla yields result a; this result controls a subsequent unitary U a , applied to the probe and the upper 
ancilla, and determines a conditional measurement M a on the upper ancilla, which has result b. The two measurement results 
then control a unitary U a ,b applied to the probe and a conditional measurement M a ,b on the probe. The left-hand circuit 
is equivalent to the circuit on the right, in which the controls are applied coherently (boxed gates Ua and Ua,b) and the 
measurements, deferred to the end of the circuit, tell one which unitary was applied. Without loss of generality, we can assume 
the measurements are described by orthogonal projectors P a , Pb\ a , an d P c \ a ,b, because any generalized measurement can be 
modeled by a projection-valued measurement on an extended system. The unitary transformations in the left-hand circuit, U a 
and U a , b, are evolution operators generated by the Hamiltonians hfho + H a (t) and h'yho + H a ,b(t), whereas the corresponding 
coherent controlled unitaries in the circuit on the right, Ua and Ua,b, are generated by the Hamiltonians hjho + Y2 a H a (t) ® Pa 
and fr"iho + ^2 a b H a ,b{t)®Pb\aPa- It is easy to verify from the evolution equations that the controlled unitaries in the right-hand 
circuit are given by Ua = Yl a ^ a ® ^ a an d Ua,b = ~}2 a b ^a,b <8> Pb\ a Pa- Thus the principle of deferred measurement can be 
rendered algebraically in the following way: if we use the left-hand circuit, the probability for obtaining results a, b, and c takes 
the form tr(C a ,b,cPoCl b ), where po is the initial state of the probe and ancillas and C a ,b,c = Pc\a,bU a ,bPb\aU a PaU; pulling 
the measurement projectors to the left in C a ,b,c changes the unitaries to the corresponding coherent controlled operations, i.e., 
C a ,b,c = P c \a,bPb\ a PaU a,bU aU , which gives the form of the probability obtained from the right-hand circuit. 



maximum (minimum) eigenvalue of H . This semi-norm 
is invariant under unitary transformations and obeys the 
triangle inequality, i.e., \\H + K\\ < \\H\\ + \\K\\ 
The importance of the semi-norm is that its square pro- 
vides a n up per bound on the variance, i.e., A 2 H < 
||F|| 2 /4 [lj. The maximum variance is achieved by pure 
states {\M H } +e" t '\m H ))/V2. 

We can now summarize the chain of inequalities satis- 
fied by the estimation precision, 

^ < y/ztfj < 2A# 7 (t) < ||Jr 7 (t)|| , (11) 

leaving us with the final task of bounding the semi-norm 
of K~{t) for the dynamics of Eqs. © and To do so, 
we define a new Hermitian operator, 



F 7 (t) = U^t)K,(t)U,(t) = iU\{t)^^- 



(12) 



which satisfies the evolution equation, dF 1 (t)/dt — 
Ui(t)hoU~,(t), with initial condition -^(O) = 
iU^{0)(dU-y{0)/d-f) = 0, since E7 7 (0) = I. Straightfor- 
ward integration provides Fy{t), and conversion back to 
Ky(t) gives 

K 7 (t)= f dsU^t)Uhs)h U 7 (s)UUt) . (13) 
Jo 

The triangle inequality and the unitary invariance of the 
semi-norm imply that 



\\KJt)\\ < / ds\\UJt)UUs)h UJs)U^t)\\ 



which gives us the desired bound on the sensitivity, 



t\\h \ 



(15) 



This bound on the estimation precision applies for any 
coupling Hamiltonian ho. It shows that the optimal sen- 
sitivity is determined by ho — indeed, it is determined by 
the range of energies in ho — and cannot be improved by 
use of a parameter-independent auxiliary Hamiltonian H 
or of ancillas not coupled directly to the parameter, al- 
though both of these might be used in physical settings to 
make the required optimal measurement accessible [14j . 

The result underlying the bound (j 1 5[) was obtained by 
Giovannetti et al. 0] for the case of discrete operations, as 
opposed to continuous time evolution, and was used there 
to show that multi-round protocols with single-system 
probes (and allowing for ancillas and adaptive measure- 
ments) have the same optimal sensitivity as single-round 
protocols with multi-system entangled probes. We use 
the bound in a different way, and our derivation shows 
directly that the ultimate sensitivity cannot be improved 
when the auxiliary Hamiltonian H acts simultaneously 
with the coupling Hamiltonian ho. 

We now apply the bound (TH)]) to draw physical conclu- 
sions about the sensitivity scaling for the various forms 
of ho- For the separable, symmetrically coupled Hamilto- 
nian of Eq. (p}, we recover the 1/N scaling of the Heisen- 
berg limit. For the symmetric /c-body coupling of Eq. f3j). 
the triangle inequality applied to the semi-norm, 



iJmll < 



< *IIM 



(14) 



E i 



\h 



(*) 



\hW\ 



(16) 
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gives a sensitivity limit that scales as 1 /N k . 

An important special case occurs when H(t) = 0, 
so that K y (t) = tho, the fc-body coupling terms in h a 

are products of single-system operators, i.e., - fc = 

hj 1 ■ ■ ■ hj k , and the single-system operators have nonncg- 
ative eigenvalues. Then the inequalities in Eqs. (fT4|) and 
(fT6|) become equalities, and an initial cat state ([2]) attains 
the maximum deviation, i.e., 

AK 7 (t) = l\\K 7 (t)\\ = i\\ho\\ = \ (X k M - Kn) ■ 

(17) 

The brief discussion of attaining the Heisenberg limit 
(k = 1), just after Eq. ©, can be applied directly to 
achieving the sensitivity limit for arbitrary k, except that 
the relative phase generalizes to cj)(t) — jt\l) [\ k M — \ k n ). 

Our result can be used to analyze the recent paper by 
Roy and Braunstein (RB) [l(J, which inspired the work 
we report here. In our notation, RB consider a system of 
iV qubits with dimensionless coupling Hamiltonian 

1 N 
h a = -(£+ + £_), Z ± = ]](X j ±iY j ), (18) 

3=1 

where Xj and Yj are Pauli operators for the j'th qubit. 
When the products are multiplied out, ho becomes a sum 
of 2 JV_1 commuting Pauli products; it has maximum de- 
viation Aha = \\h>o\\/2 = 2 which gives a quantum 
limit that scales exponentially in N. RB suggest that 
their Hamiltonian describes the atomic transitions of N 
atoms associated in a molecule, but the fundamental cou- 
pling in this case is a separable sum, as in Eq. ([1]), de- 
scribing separate transitions for each atom. The RB cou- 
pling could arise as an effective iVth-order process, but it 
would not be justified to neglect processes of other orders. 
To achieve the RB Hamiltonian as a fundamental interac- 
tion would require coupling the atoms to a rank-iV tensor 
field, but in this case, every value of N would involve a 
different fundamental coupling. One could scarcely claim 
to be estimating the same coupling constant as N changes 
if the fundamental interaction is changing. 

The most realistic possibility for taking advantage of 
multi-body couplings in parameter estimation will be for 
pairwise couplings (k = 2). Hamiltonians with sym- 
metrically parameterized two-body terms arise naturally 
in field-theoretic systems, such as quantum degenerate 
gasses, superconductors, and atomic ensembles coupled 
to a common electromagnetic field mode. Atom-atom 
interactions in a Bose-Einstein condensate (BEC) [l5| 
might offer a physically realistic approach to surpassing 
the conventional Heisenberg limit, possibly even achiev- 
ing I/TV 2 scaling. We envisage situations where an ex- 
ternal field modulates the strength of the two-body scat- 



tering term in the second-quantized condensate Hamilto- 
nian. Such a modulation occurs both for a magnetically 
tuned Feschbach resonance and for density variations due 
to gravitational gradients. 

While exponential sensitivity improvements appear 
unphysical, more modest quadratic or other polynomial 
improvements beyond the Heisenberg limit could be es- 
sential for achieving the sensitivities required in the most 
demanding precision measurements. This work was sup- 
ported in part by ONR Grant No. N00014-03-1-0426 
and AFOSR Grant No. FA9550-06-01-0178. The authors 
thank H. Barnum, A. Datta, S. Merkel, R. Raussendorf, 
and A. Shaji for helpful discussions. 
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